The P323L substitution in the SARS-CoV-2 polymerase (NSP12) confers a selective advantage during infection.

BACKGROUND: The mutational landscape of SARS-CoV-2 varies at the dominant viral genome sequence and minor genomic variant population. During the COVID-19 pandemic, an early substitution in the genome was the D614G change in the spike protein, associated with an increase in transmissibility. Genomes with D614G are accompanied by a P323L substitution in the viral polymerase (NSP12). However, P323L is not thought to be under strong selective pressure. RESULTS: Investigation of P323L/D614G substitutions in the population shows rapid emergence during the containment phase and early surge phase during the first wave. These substitutions emerge from minor genomic variants which become dominant viral genome sequence. This is investigated in vivo and in vitro using SARS-CoV-2 with P323 and D614 in the dominant genome sequence and L323 and G614 in the minor variant population. During infection, there is rapid selection of L323 into the dominant viral genome sequence but not G614. Reverse genetics is used to create two viruses (either P323 or L323) with the same genetic background. L323 shows greater abundance of viral RNA and proteins and a smaller plaque morphology than P323. CONCLUSIONS: These data suggest that P323L is an important contribution in the emergence of variants with transmission advantages. Sequence analysis of viral populations suggests it may be possible to predict the emergence of a new variant based on tracking the frequency of minor variant genomes. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

Characterisation of SARS-CoV-2 genomic variation in response to molnupiravir treatment in the AGILE Phase IIa clinical trial.

Molnupiravir is an antiviral, currently approved by the UK Medicines and Healthcare products Regulatory Agency (MHRA) for treating at-risk COVID-19 patients, that induces lethal error catastrophe in SARS-CoV-2. How this drug-induced mechanism of action might impact the emergence of resistance mutations is unclear. To investigate this, we used samples from the AGILE Candidate Specific Trial (CST)-2 (clinical trial number NCT04746183). The primary outcomes of AGILE CST-2 were to measure the drug safety and antiviral efficacy of molnupiravir in humans (180 participants randomised 1:1 with placebo). Here, we describe the pre-specified exploratory virological endpoint of CST-2, which was to determine the possible genomic changes in SARS-CoV-2 induced by molnupiravir treatment. We use high-throughput amplicon sequencing and minor variant analysis to characterise viral genomics in each participant whose longitudinal samples (days 1, 3 and 5 post-randomisation) pass the viral genomic quality criteria (n = 59 for molnupiravir and n = 65 for placebo). Over the course of treatment, no specific mutations were associated with molnupiravir treatment. We find that molnupiravir significantly increased the transition:transversion mutation ratio in SARS-CoV-2, consistent with the model of lethal error catastrophe. This study highlights the utility of examining intra-host virus populations to strengthen the prediction, and surveillance, of potential treatment-emergent adaptations.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

Analysis of SARS-CoV-2 in Nasopharyngeal Samples from Patients with COVID-19 Illustrates Population Variation and Diverse Phenotypes, Placing the Growth Properties of Variants of Concern in Context with Other Lineages.

New variants of SARS-CoV-2 are continuing to emerge and dominate the global sequence landscapes. Several variants have been labeled variants of concern (VOCs) because they may have a transmission advantage, increased risk of morbidity and/or mortality, or immune evasion upon a background of prior infection or vaccination. Placing the VOCs in context with the underlying variability of SARS-CoV-2 is essential in understanding virus evolution and selection pressures. Dominant genome sequences and the population genetics of SARS-CoV-2 in nasopharyngeal swabs from hospitalized patients were characterized. Nonsynonymous changes at a minor variant level were identified. These populations were generally preserved when isolates were amplified in cell culture. To place the Alpha, Beta, Delta, and Omicron VOCs in context, their growth was compared to clinical isolates of different lineages from earlier in the pandemic. The data indicated that the growth in cell culture of the Beta variant was more than that of the other variants in Vero E6 cells but not in hACE2-A549 cells. Looking at each time point, Beta grew more than the other VOCs in hACE2-A549 cells at 24 to 48 h postinfection. At 72 h postinfection there was no difference in the growth of any of the variants in either cell line. Overall, this work suggested that exploring the biology of SARS-CoV-2 is complicated by population dynamics and that these need to be considered with new variants. In the context of variation seen in other coronaviruses, the variants currently observed for SARS-CoV-2 are very similar in terms of their clinical spectrum of disease. IMPORTANCE SARS-CoV-2 is the causative agent of COVID-19. The virus has spread across the planet, causing a global pandemic. In common with other coronaviruses, SARS-CoV-2 genomes can become quite diverse as a consequence of replicating inside cells. This has given rise to multiple variants from the original virus that infected humans. These variants may have different properties and in the context of a widespread vaccination program may render vaccines less effective. Our research confirms the degree of genetic diversity of SARS-CoV-2 in patients. By comparing the growth of previous variants to the pattern seen with four variants of concern (VOCs) (Alpha, Beta, Delta, and Omicron), we show that, at least in cells, Beta variant growth exceeds that of Alpha, Delta, and Omicron VOCs at 24 to 48 h in both Vero E6 and hACE2-A549 cells, but by 72 h postinfection, the amount of virus is not different from that of the other VOCs.

Mutations that adapt SARS-CoV-2 to mink or ferret do not increase fitness in the human airway.

SARS-CoV-2 has a broad mammalian species tropism infecting humans, cats, dogs, and farmed mink. Since the start of the 2019 pandemic, several reverse zoonotic outbreaks of SARS-CoV-2 have occurred in mink, one of which reinfected humans and caused a cluster of infections in Denmark. Here we investigate the molecular basis of mink and ferret adaptation and demonstrate the spike mutations Y453F, F486L, and N501T all specifically adapt SARS-CoV-2 to use mustelid ACE2. Furthermore, we risk assess these mutations and conclude mink-adapted viruses are unlikely to pose an increased threat to humans, as Y453F attenuates the virus replication in human cells and all three mink adaptations have minimal antigenic impact. Finally, we show that certain SARS-CoV-2 variants emerging from circulation in humans may naturally have a greater propensity to infect mustelid hosts and therefore these species should continue to be surveyed for reverse zoonotic infections.